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Abstract 

The Cosmic Microwave Background (CMB) fluctuations at very small angular 
scales (less than 10') induced by matter sources are computed in a simplified way. 
The result corrects a previous formula appearing in the literature. The small scale 
power spectrum from cosmic strings is then calculated by a new analytic method. 
The result compares extremely well with the spectrum computed by numerical 
techniques (when the old, incorrect, formula is used). The upper bound on the 
string parameters derived from OVRO data is re-examined, taking into account 
the non-Gaussian nature of stringy perturbations on small scales. Assuming a 
conventional ionization history, the bound is jGfi < 11 x 10~ 6 , where j 2 is the 
number of horizon lengths of string per horizon volume. Current simulations give 
7 2 = 31 ± 7. 
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1 Introduction 



It is now widely accepted that cosmic microwave background (CMB) fluctuations test 
the physics of the very early universe. With the recent flurry of experiments on vari- 
ous angular scales (Meyer, Cheng & Page 1991; Smoot et al. 1992; Gaier et al. 1992; 
Meinhold et al. 1993; Ganga, Cheng, Meyer & Page 1993; Schuster et al. 1993; Gunder- 
sen et al. 1993), we are becoming more able to eliminate theories of the origin of these 
fluctuations. There are two main contenders in the competition to explain the origin of 
these fluctuations: quantum fluctuations in the metric during inflation (see for example 
Brandenberger, Feldman & Mukhanov 1992), and various kinds of semiclassical dynam- 
ics after phase transitions in field theories (Kibble 1976; Zel'dovich 1980; Vilenkin 1980; 
Turok 1989; Barriola & Vilenkin 1989; Bennett & Rhie 1990). Both are able to generate 
a more or less Harrison-Zel'dovich spectrum of density fluctuations, and both produce 
CMB fluctuations which are consistent with the COBE observations (Smoot et al. 1992), 
as far as the theoretical uncertainties allow (Peebles 1982; Bond & Efstathiou 1987; Ben- 
nett, Stebbins & Bouchet 1992; Bennett & Rhie 1993; Pen, Spergel & Turok 1993). The 
physics is quite different in each case, and each suffers in different ways from theorists' 
prejudices. However, it is to be hoped that future measurements of the spectrum of CMB 
fluctuations on a wide variety of angular scales will be able to distinguish experimentally 
between theories. 

This paper is concerned with calculating the very small angular scale fluctuations 
from cosmic strings (Bouchet, Bennett & Stebbins 1988). Strings are perhaps the most 
venerable of the theories based on the dynamics of field theories during and after phase 
transitions. Despite predating the inflationary scenario, the theory has suffered from 
analytic intractability and is consequently comparatively undeveloped. Early work on 
string seeded galaxy perturbations (Brandenberger & Turok 1986; Traschen, Branden- 
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berger & Turok 1986) was based on the first numerical simulations (Albrecht & Turok 
1989) whose detailed results have not all been confirmed by the two subsequent groups 
to work on the subject (Bennett & Bouchet 1989; Bennett & Bouchet 1990; Allen & 
Shellard 1990a; Allen & Shellard 1990b). What is lacking is a good analytic understand- 
ing of the evolution of the string network, although progress is being made in this regard 
(Kibble & Copeland 1990; Copeland, Kibble & Austin 1992; Embacher 1992a; Embacher 
1992b). It is the intention of this work to try and create an analytic approach to the 
calculation of CMB fluctuations from strings. Although the results presented apply only 
to fluctuations on very small angular scales (less than about 10'), these are precisely the 
scales on which strings make their most distinctive contribution. 

The philosophy behind the current approach is quite simple. It is that statistical 
measures of the string network itself can be used to calculate the statistics of the CMB 
fluctuations. One can in fact guess the form of simple string correlation functions on 
general grounds, backed up by some intuition gained from numerical simulations. In 
this paper, Stebbins' (1988) formula for small angle CMB anisotropies is rederived in 
a more direct way, with an error in his and previous versions of this work corrected 
(Stebbins 1993). Then, the two-point string correlators are used to calculate the very 
high frequency end of the power spectrum of CMB fluctuations. The result is compared 
to a numerical computation of the spectrum by Bouchet, Bennett and Stebbins (1988), 
using the old anisotropy formula. In view of the approximations made, it is gratifying 
that the current approach reproduces both the shape and the amplitude of the spectrum 
very well. Limits on the string linear mass density are then re-examined. The results 
can only be trusted on very small angular scales, for which the best experimental limits 
currently available are derived from recent VLA observations (Fomalont et al. 1993) 
and from OVRO (Readhead et al. 1989; Myers, Readhead& Lawrence 1993). The best 
experimental geometry for finding string is the recent RING experiment at OVRO (Myers, 
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Readhead& Lawrence 1993), which consists of 96 overlapping double difference fields in 
a circle of radius ~ 2° around the North Celestial Pole. However, sky coverage has been 
increased at the expense of sensitivity, so the best limits on the r.m.s. fluctuations still 
come from the NCP observations (Readhead et al. 1989). Other experiments at larger 
scales have been used in the past to constrain the string scenario, although in view of the 
current uncertainty surrounding the predictions of strings these should not be regarded 
as reliable. 

The plan of the paper is as follows. Firstly, I present the rederivation of Stebbins' 
(1988) formula, for small angle microwave anistropies. I then use this formula to derive an 
expression for the CMB power spectrum, subject to some plausible assumptions. These 
assumptions are justified by comparison to the numerical work of Bouchet, Bennett 
and Stebbins (1988), henceforth referred to as BBS. I then use the derived correlation 
function to obtain an upper limit on a certain combination of string parameters. To do 
this a revised Bayesian analysis is performed on OVRO data, which models the non- 
Gaussian statistics of strings. The bound, presented in ( p.!9| ), is expressed as one on 



7<j/x, where 7 2 is the number of horizon lengths of string per horizon volume, for 7 is 
not well determined as yet. Furthermore, this combination appears in all calculations of 
string-induced perturbations, and so bounding it is more useful than simply bounding 
the string tension. 

2 CMB anisotropics from moving strings 

In this section I derive the formula for the temperature pattern induced in the CMB 
by a string moving in front of it. A Minkowksi space background is used, which limits 
the applicability to fluctuations inside the horizon at decoupling. Stebbins (1993) has 
corrected his original formula (Stebbins 1988), which was also incorrectly derived in 
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earlier versions of this work. We are now, happily, in agreement. 

Imagine a large box of length b and cross-sectional area A. In this box is some 
string with spacetime coordinates X M (<r, r), where a and r are respectively spacelike and 
timelike worldsheet coordinates. Consider the gravitational effect of the string on a set of 
photons passing through this box with momentum p^ = (E; 0, 0, E). Their unperturbed 
geodesies are 

Z^x^ + Xp". (2.1) 

We would like to know what their energies are at x 11 = (to; x, z ) as a result of the string's 
gravitational field. To first order in linear theory we have 

8p,{x) = -\ [ X0 dX K^(Z(X))p^, (2.2) 

J Ai 

where h^ v is the perturbation around the Minkowski metric rj^ = diag (1,-1,— 1,-1). 
To go further it is very convenient to choose the harmonic gauge, or 

hT, v - \h» = 0, (2.3) 

where h = h u u . In this gauge the first order field equations become 

<9 2 V = 167tG(T^ - |^T), (2.4) 

where T^ u is the stress tensor. The string has its own gauge freedom, which we restrict 
by choosing the conformal gauge, or 

X-X = 0, X 2 + X 2 = 0, (2.5) 

where the dot and the prime denote differentiation with respect to r and a respectively. 
In this gauge the stress tensor is given by 

T*™ = nf drda{X^X u - X>*X U )5 (4) (x - X(a, r)) . (2.6) 
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One way to proceed from here is to compute h pv from T^ v with a retarded Green's func- 
tion, differentiate to get h up ^, and finally integrate with respect to the photon affine 
parameter (Stebbins 1988). However, a great simplification is afforded by instead calcu- 
lating 'V 2 5p p , where V is the partial derivative with respect to the transverse coordinates 
x. This is done in a roundabout fashion, using the fact that (dt — d z )f(Z) = E~ 1 df/dX: 



V 2 % = {% - % - ffl)6 Pi 



± i Xo d\ [{^ + d z )(d t - d z ) - 9 2 ](/i Wft /jy) 
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Ai 



° -SiiGdJ Xa d\T vp p v p p . (2.7) 



Ao 
Ai 

For temperature fluctuations we are interested in the variation on the photon energy 
5po = 5E. For this component we may use a cunning identity due to Stebbins (1993), 
whose proof is reproduced in the Appendix. Writing p p = p^/E, it is 

fT vpfl = -V\T ip - \-^fT ip . (2.8) 

Here, V\ = (^ — p t p')dj is the transverse derivative. Thus we have 

,SE „ r x ° 



V 2 — = 8ttG / Ed\V\T ip f 
E JXi 



+ \ 



(d t + d z )d t h - lGTrGfpPTJ X ° . (2.9) 



The terms in the square brackets are fluctuations on the bounding surfaces of our imag- 
inary box. Those at the observing time are negligible, since we are far from any sources, 
but if the initial surface is the decoupling time, there will be important fluctuations 
present. Our justification for dropping these is the finite thickness of the last scattering 
surface, which acts to smear out fluctuations on scales greater than about 10', assuming 
that z dec ~ 1000. 

We now evaluate (|2.9| ) for string sources. It is customary in the string literature 
to choose the temporal gauge X° = r; that is, identify the worldsheet time with the 



background time coordinate. In that case we find (using the fact that 5T/T = 5E/E) 

V 2 ^ = -SvrG/i jda[x- jf^j*) ' ^ 2 \x - X), (2.10) 

where worldsheet variables are evaluated at t T , given by x + = X + (cr,t T ), or t T = t + z — 
X 3 (a, t r ). Once again, the reader is cautioned that this formula, again due to Stebbins 
(1993), corrects Stebbins' (1988) result, and earlier versions of this work. 
The expression is simplified if instead we use the light cone gauge 

X + (a,r) = r. (2.11) 

The time parameter r then labels the intersections with a set of null hyperplanes with 
the worldsheet, which is in hindsight an intuitive thing to do, since our photon geodesies 
= x^ + form just such a set. Then we find 

ST r 

V 2 — = SnGfi J daX ■ V5 (2 \x - X), (2.12) 

where quantities are now evaluated at r = x + . 

These is a nice result: it says that in Minkowski space the temperature distortions 
caused by moving strings depend only on the apparent positions of the strings and their 
light cone gauge transverse peculiar velocities, and not on the entire history of the world- 
sheet. In using this gauge we will need to exercise care in taking results from numerical 
simulations, where correlation functions are always measured in the temporal gauge. 

3 Power spectrum 

In this section an expression is derived for the power spectrum of temperature fluctuations 
in terms of the two point correlators of the transverse coordinates of the string, X A (a). 
The basic two point functions are 

(X A (a)X B (a>)), (X A (a)X B (a>)), (X A (a)X B (a>)). (3.1) 
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The angle brackets denote an average over an ensemble of strings in our imaginary box. 
The starting assumption is that the string ensemble is a Gaussian process: that is, all 
correlators can be calculated in terms of the two point functions. This has not been 
tested directly, but the results can be regarded as justifying the means. 

We now make some more justifiable assumptions about the ensemble: (i) rotation, 
reflection and translation invariance of the transverse coordinates; and (ii) worldsheet 
reflection and translation invariance. Then we can reduce the number of correlation 
functions to four: 

(X A (a)X B (a>)) = \5 AB V{o-o>), (3.2) 
(X A (a)X B (a')) = ±6 AB M 1 (<T-<j') + ±e AB M 2 (<T-o'), (3.3) 
(X A (a)X B (a')) = \5 AB T(a-a r ). (3.4) 

For later convenience two other correlators will be defined: 

r(a-a') = I" do, f d<j 2 T{a 1 -a 2 ) = ({X A { ( j)-X A {<j')f), (3.5) 
U(a-a') = [* da 1 M 1 (a 1 -a) = ((X A (a)-X A (a'))X A (a)). (3.6) 

J a 1 

V and T are symmetric in their argument, while Mi and M 2 are antisymmetric. It turns 
out that we will not need the mixed correlator M 2 . 

The two point correlation function of the temperature fluctuations is defined to be 

C(r) = (ST(x)ST(x + r))/T 2 (3.7) 

where again the angle brackets denote an average over the string ensemble. In terms of 
the two dimensional Fourier transform 

6 k = j d 2 x^(x)e ik - x , (3.8) 

we have 

c ^=jlwr^ e ~ ik ' r - (3 - 9) 
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The power spectrum is the ensemble average of |5^| 2 - The Fourier transform of the 
temperature fluctuation formula ( 2.12|) is 



- k 2 5 k = i87rG^k A J daX A (a)e tk ' X ( a \ (3.10) 

Thus, normalising for the moment to unit box area A, 

P(k) = (SttG/x) 2 ^ J dado' (x A {a)X B {a')e lk ' " X (°'))} . (3.11) 

With our assumptions about the string correlation functions the ensemble average can 
be reduced to 

P{k) = |(87rG7i) 2 -^ J dado' (v(a - a') + \k 2 Tl 2 (a - a')) e -^{°-°')l±. ( 3 . 12 ) 

At this point one should check the Gaussian approximation by measuring V, Mi and T in 
one's favourite string simulation, computing P(k) from ( |3.12 ), and then comparing with 



1 5fc | 2 computed by a direct Fourier transform of the right hand side of ( ^.1U| ) averaged 
over a set of string configurations {X(cr), X(a)}. The latter procedure was adopted by 
BBS using a string simulation with the code of Bennett and Bouchet (Bennett & Bouchet 
1989; Bennett & Bouchet 1990). Unfortunately, there is very little correlation function 
data available, so we will have to content ourselves with making some educated guesses 
for the forms of V, Mi and T, based on visual inspection of the simulations. 

First, however, we must make the connection between the Minkowski space formal- 
ism and the real behaviour of photon geodesies in an expanding universe. Suppose the 
universe is flat with background metric g^ u = a 2 (7])i]^ u , where t] is the conformal time. 
Our calculations are valid for comoving momenta and comoving coordinates, and scaled 
metric h^ u = 5g^ u /a 2 , provided we can neglect terms of order h^d/a in comparison to 
hfj, v . This means that our approximation can only be trusted for perturbations on scales 
inside the causal horizon at the time of decoupling of matter and radiation, ^ dec . Bear- 
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ing that in mind, we would like to infer things about our light cone gauge correlation 
functions in a comoving box of area A and length t]q — f]dec- 

Let us outline what is believed about the three dimensional two point correlators, 

(X\<j)X\a')) = \G + {<j-<j>), (3.13) 
{X\a)X\a')) = iG-(a-a'), (3.14) 
(X\a)X\a>)) = \n(a-a'), (3.15) 

where Embacher's (Embacher 1992a; Embacher 1992b) notation has been echoed. The 
correlators are all dimensionless functions, and so must have a scale £ in them if they are 
to be non-trivial. The evidence from the numerical simulations is that this correlation 
length is proportional to the most important physical scale in the expanding universe, 
the horizon size 77. Let s = (a — a')/^(ri). Then it is found that 

i 2 f d Sl f ds 2 G + ( Sl - s 2 ) ~ £ V. (3.16) 
Jo Jo 

The exponent v is a function of scale: 

1,00- a ifs>>1 ' (3.i7) 

LI as s — > 0. 

That is, on large scales, the network behaves as a Brownian random walk, while on very 
small scales it is approximately straight. There is numerical evidence for an "intermediate 
fractal" region for s £ 1 (Bennett & Bouchet 1990; Allen & Shellard 1990b), where the 
exponent slowly interpolates between 1 and 0.5. 

The velocity correlation function G~ starts out at 2v 2 , twice the mean square velocity, 
and vanishes rapidly for s ^> 1 as befits a random walk. There is also a constraint 

/ dsG-(s) = (3.18) 
Jo 

which arises because there can be no coherent velocities in the network above the horizon 
scale. Thus G~(s) must go negative somewhere. The form of the mixed correlator Q 
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is also subject to the same integral constraint, and must also vanish at s = by its 
antisymmetry. The inferred forms of the correlators are displayed in Figure 1. 

The light cone gauge correlators contain string coordinates evaluated at all times 
between r]^ ec and 770, and thus have contributions from the equal time temporal gauge 
correlators at all times in this interval, as well as from unequal time ones. However, 
most of the string in the box is near the initial time surface, so we will assume it is the 
correlation length of the string network at r/dec that sets the scale for the light cone gauge 
correlators. 

To check that this is true, we note that the physical length density of the string is 
proprtional to £ -2 , a result of having of order one segment of length £ in a volume £ 3 . 
If we write £ =77/7, where 7 is the conformal time, then r y 2 is the number of physical 
horizon lengths of string per horizon volume rj 3 , and the comoving length of string in 
comoving volume Adrj is 



Thus most of the string is near r/dec, and it is justifiable to assume that the light cone 
gauge correlators T, V and Mi have the same general form as G + , G~ and Q respectively, 
with scale £(^dec)/a(^dec)- The exponent v is likely to be greater than 0.5, for the apparent 
position of the string need not be precisely Brownian: the nearer sections of string are 
straighter. 

We are now in a position to derive the asymptotic behaviour of P(k) as (k£) gets 
both large and small. Let us examine the first term in the power spectrum ( |3.12|) , 




(3.19) 



where rjo is toay's conformal time. In the radiation era, a(rj) = (rj/rjo) 2 , so 




(3.20) 




(3.21) 
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where o± = a ± a'. For ^ > 1 we find 

1 / Air \ 2 1 
P v (k) c (SnGrf— Lv 2 [ — - k , (3.22) 

where L is the total length of string in the box. For i( « 1 we define £ such that 
r(er) — > £|er| at large <r, and find 

P v {k) ~ ( 8 7rG/i) 2 ^(LO^iA: 2 e 2 , (3.23) 

where -Ux is a constant related to a moment of V: 

u 1 = --^J\s\V(s)ds. (3.24) 

The contribution to the power spectrum from the mixed correlator Mi is 

P M {k) = (8irGn) 2 - f rf ( 7+rfa_n 2 (a_)e- fe2r(<T - )/4 . (3.25) 
8 J 

IT 2 goes as at at small <r_, so 

P M (A;)~(87rG' /U ) 2 (LO(A;0~ 5 - (3.26) 

For small /c^, we define another constant %\ — — J \s\U 2 (s)ds/Av 2 to obtain 

p M (k) = li&irGiiYimhk 2 ?. (3.27) 

Thus the mixed correlator contributes little at the asymptotes of the spectrum, and the 
form of k 2 P{k) is now clear. It rises as k 2 for k and falls of as k~ x for high spatial 

frequencies. This general form is to be expected: at low resolution, the temperature 
pattern is uncorrelated, and so the power spectrum vanishes at low spatial frequency, 
while at high resolution we see a collection of randomly oriented edges. 

At this point it must be stressed that the function P(k) that we have been computing 
in this section does not represent the only contribution to the fluctuations in the apparent 
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temperature. The other contributions include Doppler shifting from scattering off moving 
electrons at decoupling, and Sachs- Wolf contributions from metric fluctuations in the last 
scattering surface and from decaying scalar perturbations created by the motion of the 
strings as they enter the horizon between decoupling and the present time. Furthermore, 
the calculation is wrong for k < i]^ c , which is a consequence of using the Minkowski 
space formalism. 

One could imagine fixing up this latter problem by computing the contributions to 
the power spectrum 5P(k,i]) from different times, cutting off the spectrum for k < 
i]' 1 , and then adding them up for all times since r] dcc . This approach was taken by 
Bennett, Stebbins and Bouchet (1992) to compute the fluctuations on COBE scales. 
This effectively assumes that there are no correlations between the string positions and 
velocities at one time and another, which perhaps not very safe in view of the correlations 
manifest in their numerical data (see Section 4). It also neglects the scalar perturbations 
induced by the strings as they fall within the horizon. One should not necessarily regard 
the fluctuations computed thus as a lower bound, for the temperature anisotropies caused 
by the scalar modes may well be correlated with those from the discontinuities across the 
string. A related approach in position space was adopted by Perivolaropoulos (1993). 
The correlation function was taken to be an incoherent sum of contributions made each 
time the universe doubles in size. The individual contributions were modelled by a 
cosine (which is essentially a geometrical factor arising from assuming that the strings 
are perfectly straight within a horizon distance), plus a sharp cut-off at the horizon. This 
form for the correlation function is however rather different from that derived in Section 
5. 

Doppler scattering occurs within the last scattering surface, and is a result of velocity 
perturbations to the electrons. These cannot be coherent above the horizon size at 
decoupling, which corresponds to an angular scale of about 2°. They are also expected 
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to be smoothed off on a scale corresponding to the thickness of the last scattering surface, 
which is about 5z ~ 150 for the conventional ionization history, where z& cc ~ 1100. This 
corresponds to an angular scale of around 10'. Thus we can trust our power spectrum 
only on scales below this, or ki] dcc ^> 0.1. 

4 Comparison with numerical simulations 

The BBS power spectrum was obtained by solving a Fourier space version of the temper- 
ature fluctuation formula for several null slices taken from a single run of a matter era 
string simulation (Bennett & Bouchet 1990), in which the scale factor at the end of the 
slice was twice its inital value. However, the incorrect formula was used, namely 



For the purposes of testing techniques espoused in this paper, the incorrect formula must 
be used to compare the theoretical prediction with the numerical result. BBS fitted the 
power spectrum with a function with four parameters: two for the asymptotic power 
laws, one for the position of the maximum, and one for the overall amplitude. Their fit 
is expressed in the integrated form 



where Ah is the comoving horizon size at decoupling. This corresponds to low and high 
k behaviour of k 1 - 7 and k~ 1 - 2 respectively. The exponents are not significantly different 
from the predicted values of 2 and -1. To illustrate this, the numerical spectrum has 
been compared to the function 




(4.1) 



with u — (1 — (X ■ p) 2 /(X ■ p) 2 )X (all quantities being evaluated in the temporal gauge). 




(4.2) 



F(k) = k 2 P(k)/(Gfi) 2 27i, 



(4.3) 
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with 



F(k) = ao(fc/fch)2 (4 4) 



fch being 27r/Ah- Trial and error gave = 65 ± 5 and a\ = 1.8 ± 0.2. Figure 2 shows 
k 2 P(k)/2ir for the numerical spectrum, the BBS fit, and the two parameter fit (|4.3|J4.4j) . 
The errors shown are the standard deviations of twelve null surfaces in the same simu- 
lation, taken at four different times over an expansion factor of v2, in three orthogonal 
directions. They are clearly highly correlated, which makes x 2 fitting pointless. 

Not only does the theory predict the form of the spectrum at the high and low 
frequency ends from the little information that we have about correlation functions, it 
also predicts the normalisation. The theoretical asymptote of F is, when we return the 
normalizing area A to its rightful place, 

^) = ^t(^)(t)- <-> 

where P = (X ), which we would like to compare to the numerical fit ao(k^/k). Here 
we run into a problem engendered by using the light cone gauge, for the numerical 
simulations do not directly give t, v or L. We have to convert to the temporal gauge, in 
which 

v 2 ^(u 2 (r t )}, f^y((X 2 (r r ))), (4.6) 

The mean square velocity V 2 = (X l X l ) in the matter era in the Bennett and Bouchet 
simulations is 0.37 ±0.02 (Bennett & Bouchet 1990). Since the temperature fluctuations 
depend only on the transverse components, (u 2 ) must be strictly less than \V 2 . The exact 
figure cannot be calculated with the information at hand, but we can make a crude guess 
by replacing (X 3 ) 2 and (X 3 ) 2 by their mean values, and ignoring all other correlations. 
We also face the problem of evaluating terms like {(X ■ p)~ m ), with p = 2,4. We cannot 
consistently use a Gaussian approximation here, since these expectation values diverge 
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due to the contribution at X 3 = — 1, where the string has a cusp moving towards the 
observer. However, we should not expect a Gaussian distribution for quantities like X 3 
which are constrained to be less that or equal to 1. Instead, we shall do the simplest 
thing, which is a series expansion in ((X 3 ) 2 ), so that 

((X-p)~ m ) ~ l±|m(m±l)((X 3 ) 2 ). (4.7) 

We can evaluate the required correlators with one more piece of information, the temporal 
gauge constraint ((X*) 2 ) + ((X*) 2 ) = 1. Then 

((X 3 ) 2 ) = i(l-\/ 2 ), <(X 3 ) 2 > = 1\/ 2 , (4.8) 

(± 2 ) = 2 -(l-V% (X 2 ) = lv*. (4.9) 

Thus we find 

<- 2 > = \V 2 (l - |(1 - V') + 1(1 - V 2 ) 2 (l + fv 2 )) . (4.10) 

Putting in the numerical value of V 2 , we have 

v 2 = 0.18 ±0.01, (4.11) 
t 2 = 0.42 ±0.02. (4.12) 

The total projected comoving length of string per unit area is 

Q)=A h 1 [l-(^ dC c) 1/2 )h 2 (4.13) 

where 7 2 = p s \1/ f^z dcc is the number of horizon lengths of string per horizon volume 
(Copeland, Kibble & Austin 1992), which is 31 ± 7 in the matter era simulations of 
Bennett and Bouchet (Bennett & Bouchet 1990). The BBS spectra have z/z& ec = 1/2, 
so 

L\ h /A ~ 9.1 ± 2.1. (4.14) 
15 



Thus the asymptotic theoretical spectrum is 

F th (k)~(72±17)(k h /k), (4.15) 

which is in remarkably good agreement with the numerical spectrum. The assumptions 
that have been made amount to approximating the source of the small scale fluctuations 
by uncorrelated straight segments of moving string, which we see reproduces the spectrum 
very well. The error quoted in ( [4.15 ) is statistical, and mostly due to fluctuations in the 



string density between different simulations. 

Lacking data for V(s), no prediction can be made for the low frequency end. Instead 
a hostage to fortune can be created by inferring a value for u\ through comparison of the 
low frequency end of the fitted spectrum fl4.4|) with its theoretical form (|3.23|) . We find 



Ul = V^/tiCk^Y. (4.16) 

The correlation length £ can be estimated from Fig. 2 of Bennett and Bouchet (1989) to 
be about 0.2Ah in the matter era. Thus the prediction for m is 

ui ~ 0.1, (4.17) 

whose confirmation (or otherwise) will have to await more detailed measurements of 
string correlation functions. 

5 Comparison with observation at small angular scales 

The best experiment which picks out the small scales where the theoretical spectrum can 
be trusted is OVRO (Readhead et al. 1989; Myers, Readhead& Lawrence 1993), which is 
a double difference experiment, with FWHM 108" and beam throw 7'. 15. There are also 
some recent VLA observations (Fomalont et al. 1993), which are of two regions about 
7' across with up to 10" resolution. However, the the VLA does not currently have as 
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good sensitivity as OVRO. The current theoretical uncertainties in the spectrum above 
10' make any bounds derived from larger scale experiments unreliable, and it is difficult 
to even quantify the level of unreliability. 

The correlation function of the temperature fluctuations from strings is given by 

C(r) = (GrfJ^- 2 F(ky k - r (5.1) 

This is observed by apparatus with a finite resolution. If B(x) denotes the beam response, 
then the observed temperature distribution is 

T ohs (x) =T*B = J d 2 yT(x + y)B(y), (5.2) 

and so the observed correlation function must also be convolved with the beam response: 

C ohs = C*(B*B). (5.3) 

The OVRO beam response is well approximated by a Gaussian (Readhead et al. 1989), 

B*B = — ^ exp{-r 2 /4r 2 ), (5.4) 

where r = 0'.4247. We denote the correlation function smeared on a scale r by C(r , r). 
The experiment measures the difference in temperature between a central beam (Tm) 
and the average of two flanking ones (T Ri and T R2 ), separated by r s = 7'. 15. The signal 
is therefore a temperature difference 

AT = T M -|(T R1 +T R2 ), (5.5) 

and so the the mean square temperature fluctuation measured by this type of experiment 
is 

((AT/T) 2 ) = ^C(r , 0) - 2C(r , r s ) + ^C(r , 2r„). (5.6) 
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There is a scale length ( = l/k^ai in the power spectrum, which at ~ O.lAh ~ 10' 
is about half the three dimensional correlation length of the string network, as defined 
before Eq. ( |3.23| ). Defining a dimensionless wavenumber k = k(, we find 

C(r , r)= [ ^F{k)Ukt/C) exp(- K 2 r 2 /C 2 ), (5.7) 

where Jo is the zeroth order Bessel function. We take the power spectrum to be 



We must now give the spectrum its correct normalization, using ( |2.9|) . Hence we have 

A = 16v/ri((X - (X -p)/(X ■P)Xf)[l - (z/z dec )^ 2 , (5.9) 

which we can evaluate as before, including all the string between the last scattering and 
the present day. The quantity in the angle brackets turns out to be 0.37 ± 0.02, so we 
arrive at 

A) = (16.2 ± 1.3)7 2 . (5.10) 

The form (|5.8| ) is very convenient, because an approximation to the integration, valid 
when (?"o/C) 2 ^ 1; can be found in tables. Plugging in this function into the expression 
for C(ro,r), we find 

C(r„,r) ^ C(0,r) + |C"(0, r)(r /C) 2 , (5.11) 

where 

C(0,r) = — exp(-r/C). (5.12) 

Thus 

((AT/T) 2 )(G^)- 2 ^ ~ (l - ~f rsK + ie" WC ) + 0{rl/e)- (5.13) 

For small r s /(, as is relevant in a reionized universe, the right hand side is a factor 
r s /C down on the total mean square fluctuation. One can interpret this factor as the 
probability that the beam pattern will straddle a string. 
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Putting in the experimental value for r s , the theoretical value of Aq, and the fitted 
value of di, we find that the expected r.m.s. temperature fluctuations for OVRO are, 
with z dcc = 1100, 

((AT/T) 2 ) ~ (6.3 ± 0.5)( 7 G» 2 . (5.14) 

We recall that the error is a statistical one corresponding to la fluctuations in mean 
square string velocity. When z& cc <C 1100, the formula ( [5 . 1 3| ) gives 



((AT/T) 2 ) = A (2nr s /\ h )(Gfi) 2 ~ (7.1 ± 0.6)( 1 G t x) 2 (z dcc /1100)2 . (5.15) 

The OVRO NCP 95% confidence limit on the residual sky variance is (Readhead 
et al. 1989) AT sky < 58 fiK, or A sky = AT sky /T sky < 2.1 x 10" 5 with T sky = 2.73. 
This was obtained using a Bayesian analysis with a uniform prior on A sky , assuming 
that the sky fluctuations were uncorrelated between the fields and Guassian. The first 
assumption is good for string sources too, since the NCP fields are separated by 30' 
and the scale length of the correlation function at short distances is £ ~ 10'. However, 
the non-Gaussian nature of string-induced fluctuations on small scales means that we 
must repeat the analysis, using a stringy probability distribution P s for AT. This is well 
approximated by an exponential function (Gott et al. 1990; Bennett, Bouchet & Stebbins 
1993). The theoretical analysis of Moessner, Perivolaropoulos and Brandenberger (1994) 
also seems to support this contention. Thus the probability density for measuring A in 
an observation with noise a is 

P(A,a) = ^aA^]- 1 / dA' exp(- v / 2| A - A'|/A sky + A'/2a 2 ). (5.16) 

The likelihood function for the seven uncontaminated measurements {Aj, dj}, which can 
be found in Table 4 of Readhead et al (1989), is then 

L^A^XA^) = n^(A^) (5.17) 
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This function, normalized to a maximum value of 1, is plotted in Figure 3. Assuming a 
uniform prior, as do Readhead et al (1989), it gives the Bayesian probability density for 
A s k y . From it we can infer that 



which is conventionally interpreted as a 95% confidence limit. We note that this limit 
is significantly looser than the bound on Gaussian fluctuations. This is because P s is 
strongly peaked near AT = 0. 

We can now bound 'yGfi: using our 95% confidence limit ( |5.18| ), and the central 
theoretical prediction ( |5.14D , we obtain (for a universe without reionization) 



We cannot rigorously bound the string tension, for the probability distribution of G/j 
depends crucially on how the probability distribution of 7, say P 7 , behaves near 7 = 0. 
The probabilities for A and 7 are independent, so 



If P 7 does not vanish sufficiently fast as 7 — > 0, then P(Gfi) will not converge swiftly 
enough to have a variance. Since we do not know P 7 , it is probably not worth doing more 
than estimating a bound for Gfi. The 2a statistical lower bound on 7 is 4.1, so at a kind of 
90% confidence level, we can say Gfi < 3 x 10 -6 . One is entitled to question the accuracy 
of the estimates of the correlation functions. However, the estimation method reproduced 
the numerical power spectrum to within about 10%, which corresponds to about 5% in 
the calculation of jGfi. Thus for safety we should perhaps quote 'jGfi < 12 x 10 -6 . 
A more precise bound will have to await a thorough reanalysis of the BBS simulations 



p(A sky < 2.8 x 10~ 5 ) = 0.95, 



(5.18) 



'yGfi, < 11 x 10 



6 



(5.19) 




(5.20) 
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(Bennett, Bouchet & Stebbins 1993). However, this is still very stringent, and depends 
rather weakly on the decoupling redshift: the bound is raised by a factor (HOO/zdec) 1 ^ 4 - 
We conclude by noting that the RING experiment in particular (Myers, Readhead& 
Lawrence 1993) deserves close examination, for there will be several fields containing 
strings (if the strings are there of course). This experiment does in fact report a signal, 
which if fitted to a Gaussian distribution results in an excess variance of around 100/iK. 
Although the authors do not entirely exclude contamination by point sources, it would 
be interesting to test the hypothesis that this signal is due to string. 

6 Conclusions 

In this paper a better analytic understanding of the CMB fluctuations produced by strings 
has been arrived at, at least on small angular scales (less than 10'). There are a number 
of ingredients in the success of the approach, which unfortunately make its extension 
to larger scales difficult. Principally, the Minkowski space approximation results in a 
very simple formula for the projected temperature pattern, which depends only on the 
positions and velocities of the strings on the backward light cone of the observer. This 
can only be justified for fluctuations on angular scales less than a degree. An encouraging 
success is that the Gaussian approximation for strings, in which only the two point string 
correlators are used, reproduces accurately the numerical spectrum computed by BBS. 
Using the (corrected) theoretical spectrum, a limit ( |5.19|) can be derived on ^Gfi from 



observational data on small angular scales, where 7 is the number of horizon lengths 
of string per horizon volume. Because of the corrected formula, this translates to a 
more stringent limit than that given in Bennett, Bouchet & Stebbins (1989), even when 
statistical and theoretical uncertainties are properly accounted for. We also find that the 
correlation function on very small scales is approximately C(r) ~ C(0) exp(— r/Q, with 
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£ ~ 10' (see Eq. [5 . 1 2| ) . This is rather different from the Standard Cold Dark Matter 
form C(0)/(1 + r 2 /2a 2 ), with a ~ 10' (Bond & Efstathiou 1984). In experimental 
papers, a Gaussian autocorrelation function is often used to derive limits on temperature 
fluctuations, which is a reasonable approximation to the SCDM correlation function. 
However, when limits on cosmic strings are required, an exponential correlation function 
is to be preferred. 

I am extremely grateful to Albert Stebbins for many useful discussions, and partic- 
ularly for communicating the revised formula for temperature anisotropics. I am also 
grateful to him and his collaborators David Bennett and Frangois Bouchet for making 
available their data. I thank also Paul Shellard for discussing his and Bruce Allen's 
numerical results, and Ron Horgan for a factor log e 10. This work was suported by the 
SERC. 

7 Appendix 

Our purpose is here to establish the identity 



where Z M = x M + AjA used in Section 2 to derive the anisotropy formula. This we do 
using energy-momentum conservation, for 



p"T vp>0 (Z) 



V\T ip (Z) 



1 d 
EdX 



pT ip (Z), 



(Al) 



pifli 



pi,0- 



(A2) 



Now we use the fact that when acting on functions of Z, 




(A3) 



22 



Thus 

which establishes the result, for <5 lJ — p l p> projects onto the transverse coordinates. 
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9 Figure captions 

Figure 1. The form of the three dimensional two-point string correlation functions 

G+(a) = 2(X i (a)X i (0)) (Figure la); G~{a) = 2(X i (a)X i (0)} (Figure lb); and {1(a) = 
4(X i ((x)X i (0)) (Figure lc). 

Figure 2. Numerical and analytic forms of the power spectrum of CMB fluctuations 

induced at small angular scales by cosmic strings between z dec and z dec /2. The spectrum 
is defined by F(k) = A; 2 (|5^,|) 2 /27r((j'/i) 2 , and plotted against wavenumber in units of 
fch = 27r/?7dec, where r/ dec is the horizon size at decoupling. The fitted functions are given 



in equations (O) and 



Figure 3. The likelihood function for the OVRO NCP measurements (excluding one 

contaminated field), using an exponential model for the probability distribution of string- 
induced sky fluctuations. This gives p(AT sky < 77 fiK) = 0.95, which is a bound of 
2.8 x 10~ 5 on the fractional temperature anisotropy. This is about a factor of 4/3 greater 
than the Gaussian limit. 



26 




Figure lb 




Figure lc 
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